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, The Bose-Einstein condensation (BEC) of magnetoexcitonic polaritons in two-dimensional (2D) 

electron-hole system embedded in a semiconductor microcavity in a high magnetic field B is pre- 
^ ' dieted. There are two physical realizations of 2D electron-hole system under consideration: a 

graphene layer and quantum well (QW). A 2D gas of magnetoexcitonic polaritons is considered in 
a planar harmonic potential trap. Two possible physical realizations of this trapping potential are 
assumed: inhomogeneous local stress or harmonic electric field potential applied to excitons and 
a parabolic shape of the semiconductor cavity causing the trapping of microcavity photons. The 
effective Hamiltonian of the ideal gas of cavity polaritons in a QW and graphene in a high magnetic 
field and the BEC temperature as functions of magnetic field are obtained. It is shown that the 
' effective polariton mass McB increases with magnetic field as B^''^. The BEC critical temperature 

■ Ti"' decreases as B~^^^ and increases with the spring constant of the parabolic trap. The Rabi 
splitting related to the creation of a magnetoexciton in a high magnetic field in graphene and QW 

^ I is obtained. It is shown that Rabi splitting in graphene can be controlled by the external magnetic 

field since it is proportional to B~^^'^, while in a QW the Rabi splitting does not depend on the 
magnetic field when it is strong. 
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I. INTRODUCTION 



o 

O . In the past decade, Bose coherent effects of 2D excitonic polaritons in a quantum well embedded in a semiconductor 
microcavity have been the subject of theoretical and experimental studies [H, To obtain polaritons, two mirrors 
placed opposite each other form a microcavity, and quantum wells are embedded within the cavity at the antinodes 
^ ' of the confined optical mode. The resonant exciton-photon interaction results in the Rabi splitting of the excitation 
I spectrum. Two polariton branches appear in the spectrum due to the resonant exciton-photon coupling. The lower 
\l ■ polariton (LP) branch of the spectrum has a minimum at zero momentum. The effective mass of the lower polariton is 
[ extremely small, and lies in the range 10^^ — 10^''^ of the free electron mass. These lower polaritons form a 2D weakly 
; . interacting Bose gas. The extremely light mass of these bosonic quasiparticles, which corresponds to experimentally 
00 ' achievable excitonic densities, result in a relatively high critical temperature for superfluidity, of 100 K or even higher. 
, The reason for such a high critical temperature is that the 2D thermal de Broglie wavelength is inversely proportional 
to the mass of the quasiparticle. 

While at finite temperatures there is no true BEC in any infinite untrapped 2D system, a true 2D BEC quantum 
phase transition can be obtained in the presence of a confining potential [3|, . Recently, the polaritons in a harmonic 
potential trap have been studied experimentally in a GaAs/AlAs quantum well embedded in a GaAs/AlGaAs micro- 
cavity [51]. In this trap, the exciton energy is shifted using a stress-induced band-gap. In this system, evidence for the 
' BEC of polaritons in a quantum well has been observed @. The theory of the BEC and superfluidity of excitonic 
polaritons in a quantum well without magnetic field in a parabolic trap has been developed in Ref. [7]. The Bose 
condensation of polaritons is caused by their bosonic character J 9 ■ 

While the 2D electron system was studied in quantum wells [9| in the past decade, a novel type of 2D electron 
system was experimentally obtained in graphene, which is a 2D honeycomb lattice of the carbon atoms that form 
the basic planar structure in graphite [10, 11]. Due to unusual properties of the band structure, electronic properties 
of graphene became the object of many recent experimental and theoretical studies [lO|, [llj, [l^, [iSl, [ij, ll5|, ll(| . 
Graphene is a gapless semiconductor with massless electrons and holes which have been described as Dirac-fermions 
[13] ■ The unique electronic properties in graphene in a magnetic field have been studied recently [H, [H, [13, [2l| . 



The electron-photon interaction in graphene was discussed, for example, in Ref. [22|. The energy spectrum and the 
wavefunctions of magnetoexcitons, or electron-hole pairs in a magnetic field, in graphene have been calculated in 
interesting works [23l. |24|] . 
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The spatially- indirect excitons in coupled quantum wells (CQWs), with and without a magnetic field B have been 
studied recently experimentally in Refs. [2^, [2^, [13, . The experimental and theoretical interest in these systems is 
particularly due to the possibility of the BEC and superfluidity of indirect excitons, which can manifest in the CQW 
as p ersistent electrical currents in each well and also through coherent optical properties and Josephson phenomena 
|29l [sol [31I [3^ . Since the exciton binding energies increase with rnagnetic field, 2D magnetoexcitons survive in a 
substantially wider temperature range in high magnetic fields [H, [H, [33, IH, 113, S, 111] . The BEC and superfluidity of 
spatially-indirect magnetoexcitons with spatially separated electrons and holes have been studied in graphene bilayer 
[40] and graphene superlattice [41, 42,]. The electron-hole pair condensation in the graphene-based bilayers have been 
studied in [4^0,1451,14^. However, the polaritons in graphene in high magnetic field have not yet been considered. 
The BEC and superfluidity of cavity polaritons in a QW without a trap were considered in [13, [3 ■ It is interesting to 
study a 2D system such as polaritons in graphene embedded in a microcavity from the point of view of the existence 
of the BEC within it. 

The purpose of this paper is to point out the existence of the BEC of the magnetoexcitonic polaritons in a QW and 
a graphene layer embedded in a semiconductor microcavity in a strong magnetic field and to discuss the condition of 
its realization. Since it was shown that the magnetoexcitons in a QW and graphene layer in a high magnetic field 
can be described by the same effective Hamiltonian with the different effective mass of a magnetoexciton ^41j, 14^ , 
we expect to obtain the similar expressions for the critical temperature of BEC for cavity polaritons in a QW and 
graphene with the only difference in the effective mass of a magnetoexciton. 

The paper is organized in the following way. In Sec.[lT|the spectrum of an isolated magnetoexciton with the electron 
and hole in a single graphene layern and QW is derived by applying perturbation theory with respect to the strength 
of the Coulomb electron-hole attraction. In Sec. IIIII the effective Hamiltonian of microcavity polaritons in graphene 
and QW in a high magnetic field along with a trapping potential is derived. In Sec. IIVI the Rabi splitting related to 
the creation of a magnetoexciton in graphene and QW in a high magnetic field is obtained. The temperature of BEC 
and the number of polaritons in Bose-Einstein condensate as a function of temperature, magnetic field and spring 
constant are calculated in Sec.|Vl Finally, the discussion of the results and conclusions follow in Sec. I VII 



II. AN ISOLATED MAGNETOEXCITON IN A SINGLE GRAPHENE LAYER AND QW 

When an undoped electron system in graphene in a magnetic field without an external electric field is in the ground 
state, half of the zeroth Landau level is filled with electrons, all Landau levels above the zeroth one are empty, and 
all levels below the zeroth one are filled with electrons. We suggest using the gate voltage shown in Fig. [T]to control 
the chemical potential in graphene by two ways: to shift it above the zeroth level so that it is between the zeroth 
and first Landau levels (the first case) or to shift the chemical potential below the zeroth level so that it is between 
the first negative and zeroth Landau levels (the second case). In both cases, all Landau levels below the chemical 
potential are completely filled and all Landau levels above the chemical potential are completely empty. In the first 
case, there are allowed transitions between the zeroth and the first Landau levels, while in the second case there are 
allowed transitions between the first negative and zeroth Landau levels (see the selection rules for optical transitions 
between the Landau levels in single-layer graphene [i^ and the analogous rules for the transitions between Landau 
levels in a 2D semiconductor 50]). Correspondingly, we consider magnetoexcitons formed in graphene by the electron 
on the first Landau level and the hole on the zeroth Landau level (the first case) or the electron on the zeroth Landau 
level and the hole on the Landau level —1 (the second case). Note that by appropriate gate potential we can also use 
any other neighboring Landau levels n and n + 1. 

It is obvious that magnetoexcitons formed in graphene are two-dimensional, since graphene is a two-dimensional 
structure. Below we show that for the relatively high dielectric constant of the microcavity, e 3> /(hvp) ~ 2 
{vf — V^at/{2h) is the Fermi velocity of electrons in gra phe ne, where a = 2.566 A is a lattice constant and t « 2.71 eV 
is the overlap integral between the nearest carbon atoms the magnetoexciton energy in graphene can be calculated 
by applying perturbation theory with respect to the strength of the Coulomb electron-hole attraction analogously as 
it was done in [s^ for 2D quantum wells in a high magnetic field with non-zero electron and hole masses (me 7^ 
and mfi 7^ 0). This approach allows us to obtain the spectrum of an isolated magnetoexciton with the electron 
on the Landau level 1 and the hole on the Landau level in a single graphene layer. The characteristic Coulomb 
electron-hole attraction for the single graphene layer is /{trB), where e is the dielectric constant of the environment 
around graphene, rs = ch/{eB) denotes the magnetic length of the magnetoexciton in the magnetic field B, and 
c is the speed of light. The energy difference between the first and zeroth Landau levels in graphene is Kvp/tb- For 
graphene, the perturbative approach with respect to the strength of the Coulomb electron-hole attraction is valid 
when e^/lers) ^ hvp/rs [33]. This condition can be fulfilled at all magnetic fields B if the dielectric constant of 
the surrounding media satisfies the condition /{ehvp) ^ 1. Therefore, we claim that the energy difference between 
the first and zeroth Landau levels is always greater than the characteristic Coulomb attraction between the electron 
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FIG. 1: The graphene sheet in the presence of the apphed electric E and magnetic B fields. 

and the hole in the single graphene layer at any i? if e 3> /{hvp) w 2. Thus, applying perturbation theory with 
respect to weak Coulomb electron-hole attraction in graphene embedded in the GaAs microcavity (e = 12.9) is more 
accurate than for graphene embedded in the Si02 microcavity (e = 4.5). This condition for perturbation theory in 
graphene is different from the 2D quantum well in GaAs, since in the latter case the energy deference between the 
neighboring Landau levels is hujc, where ujc = eB / (c^eh) is the cyclotron frequency, fXeh — We7Ti/i/(me + mh), and mg 
and rrih are the effective masses of the electron and the hole, correspondingly |33j. Therefore, for the quantum well in 
GaAs, the binding energy of the magnetoexciton is much smaller than the energy difference between two neighboring 
Landau levels only in the limit of high magnetic field B ^ e^cn1y^/{e^ft'), and perturbation theory with respect to 
weak electron-hole attraction can be applied only for high magnetic field. 
The operator for electron-hole Coulomb attraction is 



V{r) , 

er 



(1) 



where r = — r/i, and and r^, are vectors of an electron and a hole in a 2D plane, respectively. 

A conserved quantity for an isolated electron-hole pair in a magnetic field B is the generalized magnetoexciton 
momentum P [sj, [H, [s^l , which is given by 



-ihVe ~ ihVh + -(Ae - A,0 - -[B X (re - r,,)] . 
c c 



(2) 



The conservation of P is related to the invariance of the system upon the simultaneous translation of an electron 
and a hole along with a gauge transformation. In Eq. the cylindrical gauge for the vector potential is used: 

Ae(;i) = 1/2[B X r^(^h)]- 

The eigenfunction ipT of the Hamiltonian of the two-dimensional electron-hole pair in graphene in t he perp endicular 
magnetic field B, which is also the eigenfunction of the generalized momentum P, has the form [ssl. ISSl. [53| : 



t/jp (R, r) — exp 



l>(r-po) 



(3) 



where R = (r^ + rh)/2 and po = c[B x P]/{eB^). 

The wave function of the relative coordinate $(r) in Eq. ^ can be expressed in terms of the two-dimensional 
harmonic oscillator eigenfunctions (r). For an electron at the Landau level n+ and a hole at the level n_, the 

four-component wave functions are [23[ 



/ Si 



^-(r)-(^/2) 



5„ , .o+A'^ 



-l,|n-|-l(l-) \ 
^l>-|(l") 



(r) 



(4) 



|n+|,|n- 



where s± = sgn(n±). 
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In a high magnetic field, the magnetoexciton is formed by an electron on the Landau level 1 and a hole on the 
Landau level with the following four-component wave function: 



where $„j_„2(r) is the two-dimensional harmonic oscillator cigcnfunction given by 




.(r) = (2^) 



-l/2f,-|r?j.|/2_ 



sgn(m)^ 



j,|m| 



exp 



4r| 



(5) 



2r2 



(6) 



In Eq. ([5]), denotes Laguerre polynomials, m — ni — n2, n = min(ni,n2), and sgn(m)™ = 1 for m = 0. Note 
that we consider a magnetoexciton formed by an electron and a hole located in the same type of valley, e.g., in the 
point K (or K') of Brillouin zone. 

The magnetoexciton energies En_^_.n^{P) in graphene are functions of the generalized magnetoexciton momentum 
P, and in the first-order perturbation, are equal to 



(7) 



In Eq. ([7]), Enl,n- is the energy of the electron-hole pair when the electron is at the Landau level n+ and the hole is 
at the Landau level n_, and it is given by (23| 



= — V2[sgn(n+)/M-sgn(n-)/M 



while 



,n_ {P^ 



(8) 



(9) 



where |n+n_P) = 7/'p(R, r) is defined by Eq. 

We calculate the magnetoexciton energy using the expectation value of the electron-hole Coulomb attraction for an 
electron on the Landau level 1 and a hole on the Landau level 0. Neglecting the transitions between different Landau 
levels, the first order perturbation with respect to the weak Coulomb attraction results in the following expression for 
the energy of the magnetoexciton: 



Ei^„iP) = -{lOP 



e re 



1 P 



(10) 



Denoting the averaging by the 2D harmonic oscillator eigenfunctions ^'ni,n2(r) as (nrnP| . . . |nTOP)$, where n and 
m are defined below Eq. we get the energy of a magnetoexciton created by the electron and hole on the lowest 
Landau level: 



1 P 



Vir) 



1 P 



OOP 



V{r) 



OOP 



1 P 



V{r) 



1 P 



(11) 



Following [33|, it is easy to show that, for small magnetic momenta P <C fi/rs for the electron-hole Coulomb attraction 
p]) , each matrix element in Eq. (|11[) can be expressed in terms of the binding energy and the effective magnetic mass 
of the magnetoexciton formed by an electron and a hole in the quantum well with the 2D electrons and holes: 



{hmP\V{r)\fimP)^ = -5^ + 



(12) 



^Am ^^^"^ Mnm{B) arc the binding energy and the effective magnetic mass of the magnetoexciton, respectively, corre- 
sponding to the magnetoexciton in the state with quantum numbers h and m. 

Substituting Eq. (|12p into Eq. (jlip , we get the dispersion law of a magnetoexciton for small magnetic momenta 



1 
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Moo(S) Moi(B) 
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(13) 
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Eq. (fT3|) can be rewritten in the form: 

i?i.o(P) = -4'^ + , (14) 

where the binding energy E^'g and the effective magnetic mass ijib of a magnetoexciton in graphene with the electron 
on the Landau level 1 and the hole on the Landau level are 

_ _\ ( pib) , jj. ^ Mb) , 



1 1 



ms 2 VMoo(S) Afoi(B) 
The constants £q^{B), £q'^{B), AIoo{B), and Moi{B) depend on the magnetic field B, and are given in Ref. [sj 



(15) 



£^'^iB) 


— —£o , 


£!n\B) 




MooiB) 


= Mo , 


Mm{B) 


= -2Mo 



(16) 

where £q is the magnetoexcitonic energy and Mq is the effective magnetoexciton mass in a quantum well. These 
quantities are defined as 

— 9 

■K e 



£o w „ 

2 ers 

Mo = . (17) 

Substituting Eq. (fT6|) into Eq. p5|) gives the binding energy fg' and the effective magnetic mass of the 
magnetoexciton in a single graphene layer in a high magnetic field: 

= 4^" = 4 V 2 ^' "^^ ^ = 7^ ■ ^''^ 

We can see that the effective magnetic mass of a 2D direct magnetoexciton is 4 times higher in graphene than in a 
quantum well, while the magnetoexcitonic energy is 3/4 times lower in graphene than in a quantum well at the same 
e and B. It is interesting to mention that we obtained the effective magnetic mass of the magnetoexciton in Eq. (jlSp 
using the four-component wavefunctions of magnetoexcitons in graphene given by Eqs. ^ and ([5]). This reflects the 
specific and different properties of magnetoexcitons and, therefore, magnetopolaritons in graphene compared to the 
polaritons in a quantum well without a magnetic field [7l| . 

At small magnetic momentum (P <^ f^-frs) for measuring energies relative to the binding energy of a magnetoex- 
citon, the dispersion relation ek{P) of a magnetoexciton is quadratic: 

e,{P) = , (19) 

2mBk 

where m^fc is the effective magnetic mass that depends on B and the magnetoexcitonic quantum numbers k = 
{n+, n_} for an electron at Landau level n+ and a hole at level n_. 

It is easy to see that the results for the binding energy and effective magnetic mass of the exciton with the electron 
on the Landau level and the hole on the Landau level —1 will be exactly the same as for the exciton with the 
electron on the Landau level 1 and the hole on the Landau level 0. 

We have derived above the spectrum of the single magnetoexciton in graphene ()14|) . which is described by the 
eigenfunction of Dirac equation that has the four-component spinor structure given by Eq. ([5]). Alternatively, the 
wave function of the magnetoexciton in a QW has the one-component structure, because this wave function is the 
eigenfunction of Schrodinger equation. However, Eq. (jM]) is valid also for a QW, but the binding energy and effective 
magnetic mass of 2D magnetoexciton formed by the electron and hole in the QW on the zeroth Landau level are given 

by Bi 

4'^-^^o-^/| — , mB = Mo^^^^. (20) 

V 2 ers \/i^e/rB 

Also for a QW the expression for the single magnetoexciton spectrum given by Eq. (jl9p is valid. 
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III. THE EFFECTIVE HAMILTONIAN OF TRAPPED MICROCAVITY POLARITONS IN 
GRAPHENE AND IN A QW IN A HIGH MAGNETIC FIELD 

Polaritons are linear superpositions of excitons and photons. In high magnetic fields, when magnetoexcitons may 
exist, the polaritons become linear superpositions of magnetoexcitons and photons. Let us define the superposi- 
tions of magnetoexcitons and photons as magnetopolaritons. It is obvious that magnetopolaritons in graphene are 
two-dimensional, since graphene is a two-dimensional structure. The effective Hamiltonian of magnetopolaritons in 
graphene and a QW in the strong magnetic field is given by 

P^tot — ^mex -^ph ^~ ^rnex—ph ; (21) 

where Hph is a photonic Hamiltonian, and He^c-ph is the Hamiltonian of magnetoexciton-photon interaction, and 
Hmex is a effective magnetoexcitonic Hamiltonian. Let us analyze each term of the Hamiltonian for magnetopolaritons 
pTjl . It was shown in Rcf. [4l|,|42| that 2D magnetoexcitons in graphene and a QW in a high magnetic field can be 
described by the same effective Hamiltonian Hmex ■ The effective Hamiltonian of 2D non-interacting magnetoexcitons 
in the infinite homogeneous system in a high magnetic field is given by [Zll 

Hmex — £me.x{P)bpbp , (22) 

p 

where bp and bp are magnetoexcitonic creation and annihilation operators obeying the Bose commutation relations. 
For Hamiltonian ([22|) , the energy dispersion of a single magnetoexciton in a graphene layer is given by 

Smex {P) = Etand ^ 4'^ + £0 {P) ■ (23) 

Eband = e[^q = y/2hvF /rs is the band gap energy, which is the difference between the Landau levels 1 and in 

graphene defined by Eq. ([8]). is the binding energy of a 2D magnetoexciton with the electron in the Landau level 
1 and the hole on the Landau level in a single graphene layer, and £o(^') = P^/(2mB), where ms is the effective 
magnetic mass of a 2D magnetoexciton with the electron on the Landau level 1 and hole on the Landau level in a 
single graphene layer given by Eq. ^9]) . 

It can be shown that the interaction between two direct 2D magnetoexcitons in graphene with the electron on the 
Landau level 1 and the hole on the Landau level can be neglected in a strong magnetic field, in analogy to what is 
described in Ref. (33| for 2D magnetoexcitons in a quantum well. The dipole moment of each exciton in a magnetic field 
is di^2 = epo = [B X Pi, 2] /B [33|, where Pi and P2 are the magnetic momenta of each exciton and Pi, P2 ^ l/^'s- 
The magnetoexcitons are located at a distance tb from each other. The corresponding contribution to the energy 
of their dipole-dipole interaction is - (rs/Rf PiP2r%/e ~ {rs/ Pf PiP2/{eMo) < eV|/(eP3). Inputting the 
radius of the magnetoexciton in graphene rp^i ~ tb [40| . we obtain that the van der Waals attraction of the exciton 
at zero momenta is proportional to ~ (ro,i/P)^ ~ {rs/R)^- Therefore, in the limit of a strong magnetic field for 
a dilute system tb ^ R, both the dipole-dipole interaction and the van der Waals attraction vanish, and the 2D 
magnetoexcitons in graphene form an ideal Bose gas analogously to the 2D magnetoexcitons in a quantum well given 
in Ref. [3^. Thus, the Hamiltonian pT|) does not include the term corresponding to the interaction between two direct 
magnetoexcitons in a single graphene layer. So in high magnetic field there is the BEC of the ideal magnetoexcitonic 
gas in graphene. 

Let us analyze the other two terms in the Hamiltonian (j2ip . The Hamiltonian of non- interacting photons in a 
semiconductor microcavity is given by [ssj : 

^p/. = ^ep/i(P)4ap' (24) 
p 

where cip and dp are photonic creation and annihilation Bose operators. The cavity photon spectrum is given by 



EphiP) = {c/n)^P^ + h^TT^L^^ . (25) 

In Eq. , Lc is the length of the cavity, n = y^ec is the effective refractive index and ec is the dielectric constant 
of the cavity. We assume that the length of the microcavity has the following form: 

n \rjband — <^B 
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FIG. 2: The length of the microcavity of GaAs (ec ~ 12.9), corresponding to magnetoexciton-photon resonance, as a function 
of the magnetic field B. 



corresponding to the resonance of the photonic and magnetoexcitonic branches at P = (i.e. SmexiO) = ep/i(0)). 
The length of the microcavity, corresponding to a magnetoexciton-photon resonance, decreases with the increment of 
the magnetic field as B~^/^. The dependence of the length of the microcavity corresponding to the magnetoexciton- 
photon resonance on the magnetic field is shown in Fig. [2l The resonance between magnetoexcitons and cavity 
photonic modes can be achieved either by controlling the spectrum of magnetoexcitons Sex (P) by changing magnetic 
field B or by choosing the appropriate length of the microcavity Lc- Let us mention that, while in the presence of a 
high magnetic field, the length of the microcavity corresponding to the magnetoexciton-photon resonance depends on 
the magnetic field as it is shown in Fig. [21 This effect does not take place in the system without a magnetic field 
The Hamiltonian of the harmonic magnetoexciton-photon coupling has the form [s^ : 

Hmex-ph = fi^B. ^ apbp + h.c, (27) 
p 

where the magnetoexciton-photon coupling energy represented by the Rabi constant hrtu is obtained in Sec. IIVI Let 
us mention that is obtained for a QW from the standard procedure describing the electron-photon interaction 
in the Hamiltonian by the P • A term, while in a single graphene layer flji is obtained from the electron-photon 
interaction based on the Dirac Hamiltonian for the electron in graphene. 
The excitonic and photonic operators are defined as [131 

bp — Xpp-p — Cpu-p, dp — Cppp + Xpu-p, (28) 

where pp and up are lower and upper magnetopolariton Bose operators, respectively. Xp and Cp are given by 

Xp = ^ Cp = ^ == , (29) 



2 



and the energy spectra of the lower/upper magnetopolaritons are 

£ph{P) + [P) 



£lp/up{P) = 



l^Ji£phiP) - e^exiPW + MfiW ■ (30) 



Eq. (|30p implies a splitting of 2hD,p between the upper and lower states of polaritons at P = 0, which is known as 
the Rabi splitting. Let us also mention that l-'^pp and |Cpp = 1 — |^pp represent the magnetoexciton and cavity 
photon fractions in the lower magnetopolariton. 



8 



Substituting Eq. ([28]) into Eqs. ([22]), ([241) and ^7}, we conclude that the total Hamiltonian Htot UM can be 
diagonalized by applying unitary transformations (|28D and has the form: 

£up{P)''J-pUp, (31) 

p p 

where Pp, pp, Up, up arc the Bose creation and annihilation operators for the lower and upper magnetopolaritons, 
respectively. 

Eq. ([?T|) is the Hamiltonian of magnetopolaritons in a single graphene layer in a high magnetic field. Our particular 
interest is the lower energy magnetopolaritons which produce the BEC. The lower palaritons have the lowest energy 
within a single graphene layer. Therefore, from Eq. (|3ip we can obtain 

iJtot = ^eLp(i^)p;pp . (32) 
p 

Similarly to the case of Bose atoms in a trap (55l.[56| in the case of a slowly varying external potential, we can make 
the quasiclassical approximation, assuming that the effective magnetoexciton mass does not depend on a characteristic 
size I of the trap and it is a constant within the trap. This quasiclassical approximation is valid if P ^ h/l. The 
harmonic trap is formed by the two-dimensional planar potential in the plane of graphene. The potential trap can be 
produced in two different ways. In case 1, the potential trap can be produced by applying an external inhomogeneous 
electric field or inhomogeneous local stress. The spatial dependence of the external field potential V{r) is caused 
by shifting of magnetoexciton energy by applying an external inhomogeneous electric field or inhomogeneous local 
stress. The photonic states in the cavity are assumed to be unaffected by this electric field or stress. In this case 
the band energy E^and is replaced by Efjand{i') = EbandiO) + V{r). Near the minimum of the magnetoexciton energy, 
V{r) can be approximated by the planar harmonic potential where 7 is the spring constant. Note that a high 

magnetic field does not change the trapping potential in the effective Hamiltonian [50|, l57|- In case 2, the trapping of 
magnetopolaritons is caused by the inhomogeneous shape of the cavity when the length of the cavity is given by 

= ^) V ' (33) 

n i^Eband - + iry2j 

where r is the distance between the photon and the center of the trap. In case 2, the 7 in Eq. (|33|) is the curvature 
characterizing the shape of the cavity. In case 1, for the slowly changing confining potential V{r) = the 
magnetoexciton spectrum is given in the effective mass approximation as 

e^L(^) = erne^iP) + V{r) ^ icMhTTL^' + lr^ + ^^ (34) 

where r is now the distance between the center of mass of the magnetoexciton and the center of the trap. The 
Hamiltonian for photons in this case is given by Eq. (I24p . the spectrum of photons is shown by Eq. ([25| and the 
length of the microcavity is given by Eq. (|26|) . 

In case 2, for the slowly changing shape of the length of cavity given by Eq. ([55)) . the photonic spectrum is given 
in the effective mass approximation as 



e'S(P) = + ^ [Eband ~ 4'^ + ^) . (35) 

1/2 

This quasiclassical approximation is valid if P ^ h/l, where / — (^./(tobcjo)) is the size of the magnetoexciton 
cloud in an ideal magnetoexciton gas and luq = y^'y/mB- The Hamiltonian and spectrum of magnetoexcitons in this 
case are given by Eq. ([^ and ([^ . correspondingly. 

The total Hamiltonian Htot can be diagonalized by applying unitary transformations. At small momenta a = 

l/2(m^^ + {c/n)Lc/hTT)P^ /\hQii\ ^ 1 ( Lc — hnc/n {^Eband — ^b''^ ) weak confinement (3 = ^r'^ /\h^R\ ^ 1, 
the single-particle lower magnetopolariton spectrum obtained through the substitution of Eq. (I34p into Eq. ([30l) , in 
linear order with respect to the small parameters a and (3, is 

so(P) « ^ft^P- - \m^\ + y + \ + ^) (36) 
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FIG. 3: The ratio of the effective magnetopolariton mass Mcft(i3) to the mass of a free electron me as a function of magnetic 
field B. 



Let us emphasize that the spectrum of non- interacting magnetopolaritons Eq (P) at small momenta and weak confine- 
ment is given by Eq. (j36p for both physical realizations of confinement: case 1 and case 2. By substituting Eq. (j34p 
into Eq. (|^^ . we obtain Xp w l/\/2. The condition for the validity of the quasiclassical approach in Eq. (P^ . PI ^ h, 
is also applied here. 

If we measure the energy relative to the P — lower magnetopolariton energy {c/n)hnL^^ — \hQji\, wc obtain the 
resulting effective Hamiltonian for trapped magnetopolaritons in graphene in a magnetic field. At small momenta 



a ^ 1 {Lc = hTTc/n (^Ei 



band 



-£ 



{b) 



) and weak confinement /3 <C 1, this effective Hamiltonian is 



H, 



cff 



(37) 



where the sum over P is carried out only over P 3> h/l (only in this case the quasiclassical approach used in Eq. 
is valid) , and the effective magnetic mass of a magnetopolariton is given by 

cLc{B) 



Meff(S) = 2 



(38) 



According to Eq. (|38p , the effective magnetopolariton mass M^s increases with the increment of the magnetic field as 
B^^^, as shown in Fig. [31 Let us emphasize that the resulting effective Hamiltonian for magnetopolaritons in graphene 
in a magnetic field for the parabolic trap is given by Eq. (|37p for both physical realizations of confinement represented 
by case 1 and case 2. The effective magnetic mass of the magnetopolariton in a QW is approximately the same as in 
graphene, since the contribution to AIcs{B) from the second term in the r.h.s. of Eq. ([38]) is much higher than from 
the first term. So the effective mass of the magnetopolariton in a QW can also be presented by Fig. [31 

Let us mention that the effective Hamiltonian of magnetopolaritons in a QW in microcavity is also given by Eq. (|37l) 
with the effective magnetic mass of magnetoexciton with the electron and hole on the zeroth Landau level provided 
by Eq. 



IV. THE RABI SPLITTING CONSTANT IN GRAPHENE AND A QW IN HIGH MAGNETIC FIELD 

Neglecting anharmonic terms for the magnetoexciton-photon coupling, the Rabi splitting constant Qn can be 
estimated quasiclassically as 



\hn, 



f H„ 



(39) 
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where Hint is the Hamiltonian of the electron-photon interaction. For graphene this interaction is determined by 
Dirac electron Hamiltonian as 



Hi, 



vpe 



A. 



phO 



vpe::' - 

— — a ■ £jpho , 

lUJ 



(40) 



where a = {ax,(Ty), and ay are Pauli matrices, Ap^o is the vector potential corresponding to a single cavity 

1 /2 

photon, and Epho — {8Trhxjj/{eW)) is the magnitude of electric field corresponding to a single cavity photon of 
the frequency uj in the volume of microcavity W, while for the QW this interaction is 



where 



Hint — d 12 ■ Epho , 



12 = 1 



(41) 



(42) 



is the dipole momentum of transition and the sum is taken over the coordinate vectors related to the positions of all 
the electrons in the system. 

In Eq. ([55]) the initial |i) and final |/) electron states are different for graphene and a quantum well. For the case 
of graphene these electron states are defined as 



= n<fcio)oio)i , 



l/> = hi 



(43) 



In Eq. 



is the Fermi creation operator of the electron with the y component of the wavevector k on the 



Landau level n, |0)„ denotes the wavefunction of the vacuum on the Landau level n, Y{k k\^)o corresponds to the 
completely filled zeroth Landau level, SJ^ is the Bose creation operator of the magnetoexciton with the electron 
on the Landau level n and the hole on the Landau level n' . We consider magnetoexcitons with magnetic momenta 
equal to zero, for which the Bose condensate in the system of non-interacting particles is the exact solution of the 
problem [s^l • Following Ref. [s^l Vn n' ^^^^ '^^^^ defined as 



1 



IN,, 



■ 2^ '^n',k^n,~k ' 



(44) 



where hl^, ^ is the Fermi creation operator of the hole with the y component of the wavevector k on the Landau level 

n', Nd — S/{2Trrg) is the macroscopic degeneracy of Landau levels, and S is the area of the system. 

Let us use the Landau gauge for the wavefunction of the single electron ipn,k{x,y) with the y component of the 
wavevector k on the Landau level n. In the Landau gauge with the vector potential A = (0, Bx^ 0), the two-component 
eigenf unction ipn,k{Tc) is given by [581] 



i'7i,k(x,y) 



exp(iky) 



s{n)i"- ^$„_i(a; - r|fc) 
i"<S>n{x - r%k) 



where s(n) defined by 



.in,- / (^ = 0), 



Ly are normalization lengths in the y direction, 



Cn 



1 (n = 0) , 
1/V2 (n > 0) , 



and 



$„(a;) = (2"n!y^rB) exp 



1 / X 
2 



H„ 



(45) 



(46) 



(47) 



(48) 
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where Hn{x) is the Hermite polynomiaL The corresponding eigenenergies depend on the quantum number n only 
and are given by 



2n 



(49) 



Substituting Eqs. (gl]) and (gS]) into 
from Eq. ([39]): 



evp 



and using the electron-photon interaction Hint (|40p . we finally obtain 

evp \Epho\ 



- — dx dyx 'il;l k{x,y)a ■ Epho^o,k{x,y) 



^/2( 



In Eq. ((50)) the energy of photon absorbed at the creation of the magnetoexciton (at <C £i ~ £o) is given by 



(50) 



huj ~ El — Eq — \Pi ^ 



(51) 



Substituting the photon energy from Eq. (|5ip into Eq. ()50p , we obtain the Rabi splitting corresponding to the creation 
of a magnetoexciton with the electron on the Landau level 1 and the hole on the Landau level in graphene: 



h£lvt — 2e 



Trhvprs 
V2eW 



1/2 



(52) 



As follows from Eq. the Rabi splitting in graphene is related to the creation of the magnetoexciton, which 

decreases when the magnetic field increases and is proportional to B~^^^. Therefore, the Rabi splitting in graphene 
can be controlled by the external magnetic field. Note that in a semiconductor quantum well contrary to graphene 
the Rabi splitting does not depend on the magnetic field. 

Substituting Eq. (|¥T|l and the initial |i) and final |/) electron states from Ref. [33| into ([55)1 after the integration 
we obtain the Rabi splitting constant flu for a quantum well 



where di2 is the matrix term of a magnetoexciton generation transition in a QW represented as 



/ 



(53) 



(54) 



The similar calculations for the transition dipole moment and the photon energy corresponding to the formation of 
magnetoexciton with the electron and hole on zeroth Landau level in the QW gives: 



, _ erB 

hiU = Si — Eq = huJc 



heB 



(55) 



Substituting the transition dipole moment and the photon energy given by Eq. (|55)) into Eq. we obtain the Rabi 
splitting for QW: 



Mlpi = 2eh 



1/2 



(56) 



Thus, as it follows from Eq. (|56)) . the Rabi splitting in a QW does not depend on the magnetic field in the limit of 
high magnetic field. Therefore, only in graphene can the Rabi splitting be controlled by the external magnetic field 
in the limit of high magnetic field. 

It is easy to show that the Rabi splitting related to the creation of the magnetoexciton, the electron on the Landau 
level and the hole on the Landau level —1 will be exactly the same as for the magnetoexciton with the electron 
on the Landau level 1 and the hole on the Landau level 0. Let us mention that dipole optical transitions from the 
Landau level —1 to the Landau level 0, as well as from the Landau level to the Landau level 1, are allowed by the 
selection rules for optical transitions in single-layer graphene ^4&] . 
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V. BOSE-EINSTEIN CONDENSATION OF TRAPPED MICROCAVITY MAGNETOPOLARITONS IN 

GRAPHENE AND QW 

Although Bose-Einstein condensation cannot take place in a 2D homogeneous ideal gas at non-zero temperature, as 
discussed in Ref. in a harmonic trap the BEC can occur in two dimensions below a critical temperature . Below 
we estimate this temperature. In a harmonic trap at a temperature T below a critical temperature T° {T < T^), the 
number No{T, B) of non- interacting magnetopolaritons in the condensate is given by [3] 

r(2)C(2) (.9(^)5^^) +i,i'^)g(^))Meff(i?) 
No{T,B) = N ^ 1 {kBTf 

= N- 3^2^ ^ {ksTf , (57) 

where N is the total number of magnetopolaritons, gi'^^''^'^^ and gi'^^''"'^^ are the spin and graphene valley degeneracies 
for an electron and a hole, respectively, fc^ is the Boltzmann constant, r(a::) is the gamma function and (^{x) is the 
Riemann zeta function. 

Applying the condition TVo = to Eq. ((57)) . and assuming that the magnetopolariton effective mass is given by 

Eq. (p8|) . we obtain the BEC critical temperature ri°^ for the ideal gas of magnetopolaritons in a single graphene 
layer in a magnetic field: 



^'^^"^^^^ " fc^ f (e) (e) ^ Z {h)\ .pj • (^^^ 

Igl 'gl + gl ' gl ' j Mos{B) ^ 



At temperatures above Tc^\ the BEC of magnetopolaritons in a single graphene layer does not exist. 

A three-dimensional plot of Tc^^ /^/N as a function of magnetic field B and spring constant 7 is presented in Fig.H) 
In our calculations, we used gi''^ — gi'^^ = gi'^^ = g'i^^ — 2. The functional relations between the spring constant 7 
and the magnetic field B corresponding to different constant values oiTc^^ /\/N are presented in Fig.[5l According to 
Eq. (j58p . the BEC critical temperature T^^'' decreases with the magnetic field as i?^^/^ and increases with the spring 
constant as 7^/^. These functional relations are illustrated in Figs.[4l[5]and[6l 

Substituting Eq. into Eq. ([57]) . we obtain 



13 




2 4 6 8 10 12 14 16 18 20 
Magnetic field, 6 (T) 



FIG. 5: The functional relations between the spring constant 7 and the magnetic field B corresponding to the different constant 
values of ri°VViV. We assume that the environment around graphene is GaAs with e — 12.9. 

Note that, since the quadratic spectrum of non-interacting magnetopolaritons given by Eq. does not satisfy the 
Landau criterion of superfluidity [59, 60], the ideal Bose gas of magnetopolaritons in high magnetic field in graphene 
is not a superfluid. 

Since magnetopolaritons in a QW are described by the same effective Hamiltonian as in graphene, but with the 
different magnetic mass of the magnetoexciton, the results of the calculations presented in Figs. HI [5] and [6] for the 
critical temperature of the BEC for magnetopolaritons in graphene are valid for the BEC in a QW in a high magnetic 
field. This is true because the contribution to the effective mass of the magnetopolariton from the second term in the 
r.h.s. of Eq. ([55)1 is much higher than from the first term. 

VI. DISCUSSION AND CONCLUSIONS 

In our calculations, we have assumed that the system under consideration is in thermal equilibrium. This assumption 
is valid if the relaxation time is less than the quasiparticle lifetime. Although the magnetopolariton lifetime is short, 
thermal equilibrium can be achieved within the regime of a strong pump. Porras et al. |61l| claimed that the time 
scale for polariton-exciton scattering can be small enough to satisfy this condition for the existence of a thermalized 
distribution of polaritons in the lowest fc-states in a quantum well. We expect a similar characteristic time for 
magnetopolariton-magnetoexciton scattering in graphene. However, the consideration of pump and decay in a steady 
state may lead to results which are different from the ones presented in this paper. The consideration of the influence 
of decay on the BEC may be the subject of further studies of a trapped gas. 

Above we discussed the BEC of the magnetopolaritons in a single graphene layer placed within a strong magnetic 
field. What would happen in a multilayer graphene system in a high magnetic field? Let us mention that the 
magnetopolaritons formed by the microcavity photons and the indirect excitons with the spatially separated electrons 
and holes in different parallel graphene layers embedded in a semiconductor microcavity can exist only at very low 
temperatures fc^T <C hflji. For the case of the spatially separated electrons and holes, the Rabi splitting ilji is very 
small in comparison to the case of electrons and holes placed in a single graphene layer. This is because flu ~ di2 
and the matrix element of magnetoexciton generation transition di2 is proportional to the overlapping integral of the 
electron and hole wavefunctions, which is very small if the electrons and holes are placed in different graphene layers. 
Therefore, we cannot predict the effect of relatively high BEC critical temperature for the electrons and holes placed 
in different graphene layers. 

Spin polarization is important not only for the excitations but for the condensate itself. It was shown in [47l , l48j 
that taking into account the spin degree of freedom can qualitatively modify the results for exciton-magnetopolariton 
condensation at magnetic flelds lower than the critical magnetic field. We assume that magnetic field B under 
consideration is above the critical one and, therefore, the Zeeman splitting does not affect the spectrum of collective 
excitations according to Fig. 1 in So we neglect the Zeeman splitting in our calculations. 
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FIG. 6: The ratio of the BEG critical temperature to the square root of the total number of magnetopolaritons Tc^^ as 
a function of magnetic field B at different spring constants 7. We assume the environment around graphene is GaAs with 
e = 12.9. 



To conclude, we have derived the effective Hamiltonian of the ideal gas of trapped cavity magnetopolaritons in 
a single graphene layer and a QW in a high magnetic field. The resonance between magnetoexcitons and cavity 
photonic modes can be achieved either by controlling the spectrum of magnetoexcitons Eex (P) , by changing magnetic 
field B or by controlling the length of the microcavity Lc- We analyzed two possible physical realizations of the 
trapping potential: inhomogeneous local stress or a harmonic electric field potential coupled to magnetoexcitons and 
a parabolic shape of the semiconductor cavity causing the trapping of microcavity photons. We conclude that both 
realizations of confinement result in the same effective Hamiltonian. It is shown that the effective magnetopolariton 
mass Moff increases with the magnetic field as B^^^. Meanwhile, the BEC critical temperature Ti°^ decreases as 
and increases with the spring constant as 7^/^. The gas of magnetopolaritons in graphene and a QW in a high 
magnetic field can be treated as an ideal Bose gas since magnetoexciton-magnetoexciton interaction vanishes in the 
limit of a high magnetic field and a relatively high dielectric constant of the cavity e ^ 2, according to Sec. |TT1 Let 
us mention that this condition for the high dielectric constant of the microcavity is valid only for graphene, and it 
is not valid for the quantum well. Observation of trapped cavity magnetopolaritons in graphene in a high magnetic 
field would be an interesting confirmation of the magnetopolaritonic BEC that we have described. Besides, we have 
obtained the Rabi splitting related to the creation of a magnetoexciton in a high magnetic field in graphene. Since 
this Rabi splitting is proportional to B~^/^, we conclude that the Rabi splitting in graphene can be controlled by the 
external magnetic field B, while in a quantum well the Rabi splitting does not depend on the magnetic field when it 
is strong. The results for the critical BEC temperature of magnetopolaritons in a QW and graphene in high magnetic 
field are similar, since the magnetoexcitons in both systems are described by the same Hamiltonian. 
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